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The idea of topological quantum computation (TQC) is to store and manipulate quantum information in an 
intrinsically fault-tolerant manner by utilizing the physics of topologically ordered phases of matter. Currently, 
one of the most promising platforms for a topological qubit is in terms of Majorana fermion zero modes (MZMs) 
in spin-orbit coupled superconducting nanowires. However, the topologically robust operations that are possible 
with MZMs can be efficiently simulated on a classical computer and are therefore not sufficient for realizing a 
universal gate set for TQC. Here, we show that an array of coupled semiconductor-superconductor nanowires 
with MZM edge states can be used to realize a more sophisticated type of non-Abelian defect: a genon in an Ising 
X Ising topological state. This leads to a possible implementation of the missing topologically protected tt/S 
phase gate and thus universal TQC based on semiconductor-superconductor nanowire technology. We provide 
detailed numerical estimates of the relevant energy scales, which we show to lie within accessible ranges. 


Introduction 

The promise of topological quantum computation (TQC) is 
to encode and manipulate quantum information using topo¬ 
logical qubits ma. The quantum states of a topological 
qubit do not couple to any local operators, forming a non¬ 
local Hilbert space, and are therefore intrinsically robust to 
decoherence. Currently, one of the most promising avenues 
towards developing a topological qubit is in terms of Majo¬ 
rana fermion zero modes (MZMs) in spin-orbit coupled su¬ 
perconducting nanowires 14111^. 

Topologically protected qubits are most useful if the infor¬ 
mation stored in them can be controlled in a topologically 
protected way. A system of topological qubits that are in a 
rich enough topological phase to allow complete manipula¬ 
tion of the state space would lead to the construction of a uni¬ 
versal topological quantum computer. Unfortunately, most of 
the topological phases that appear within experimental reach 
such as topological superconductors hosting MZMs lO, the 
Pfaffian fractional quantum Hall (FQH) states na, surface 
codes (TtHSI and even the recently proposed parafermion 
zero modes 126^311 are not complex enough to span the en¬ 
tire topological state space in a topologically protected way. 
In fact, the topologically protected unitary operations that are 
possible with MZMs correspond to the Clifford group, which 
can be efficiently simulated on a classical computer 1^ . Con¬ 
sequently, proposals for utilizing MZMs in quantum computa¬ 
tion require non-topological operations, and perhaps interfac¬ 
ing them with conventional, non-topological qubits 133^371 . 

Despite much previous work, it remains a major open 
problem to find a viable path towards universal TQC. Ref. 
[38l showed that some 2D systems harboring MZMs could 
support universal TQC given the possibility of dynamical 
topology changes. However, proposals to exploit these 
ideas using topological superconductors in semiconductor- 
superconductor heterostructures 1391 . or the Moore-Read 
Pfaffian FQH state with tilted interferometry l4Ql . were since 


found to be insufficient even in principle, as they lack a cru¬ 
cial ingredient: the existence of a finite energy quasiparti¬ 
cle excitation with an appropriate value of topological spin 
El. Other previously proposed physical platforms require 
either (1) exotic non-Abelian fractional quantum Hall (FQH) 
states, such as those which possess Fibonacci quasiparticles, 
or the Z 4 Read-Rezayi FQH state l^ l42U47l . or (2) complex 
designer Hamiltonians based on Josephson junction arrays, 
which effectively realize non-Abelian discrete gauge theories 
in i4^H50l . However, the viability of these proposals, even in 
principle, is questionable due to the difficulty of establishing 
the existence of non-Abelian FQH states and the impractically 
low energy scales of previously proposed designer Hamiltoni¬ 
ans El Moreover, the scientific and technological advances 
required to realize such proposals are not directly relevant to 
those being developed in the pursuit of MZMs in spin-orbit 
coupled superconducting nano wires. 

In this paper, we show that a network of coupled super¬ 
conducting nanowires hosting MZMs can be used to real¬ 
ize a more powerful type of non-Abelian defect: a genon 
||26l[3TJ|5Tl in an Ising x Ising topological state. The braiding 
of such genons can be shown to mathematically map onto the 
required dynamical topology changes of Ref. [38] and there¬ 
fore provides the missing topological single-qubit 7r/8 phase 
gate 1311 . Combined with joint fermion parity measurements 
of MZMs, these operations provide a way to realize universal 
TQC El. 

Our proposal consists of the following basic building 
blocks. First, we show that an array of suitably coupled 
MZMs in nanowire systems can realize a two-dimensional 
phase of matter with Ising topological order ca. We do this 
by showing how to engineer an effective Kitaev honeycomb 
spin model in a realistic physical system of coupled Majorana 
nano wires, where each effective spin degree of freedom cor¬ 
responds to a pair of Majorana nano wires. We present two 
approaches to doing this, corresponding to whether capacitive 
charging energies or Josephson couplings are the dominant 
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FIG. 1: (A) A pair of semiconductor-superconductor nanowire 

structures, labelled as A and B. The blue regions represent the su¬ 
perconductor, the black line represents the semiconducting nanowire 
in the topological superconducting phase, and the red dots represent 
the MZM edge states. (B) Circuit model for (a). (C) A plaquette 
of a brick lattice consisting of six effective sites. The semiconduc¬ 
tor nanowires fully extend in the horizontal direction; the gray re¬ 
gions indicate normal (non-topological) regions of the semiconduc¬ 
tor nano wires. The required applied magnetic field can be taken to be 
normal to the plane. Ia and t b parameterize the electron tunneling 
as shown. Dashed rectangle indicates a unit cell of the brick lattice, 
each containing four superconducting islands, or two pairs of A and 
B -type structures, labelled 1, • • • , 4 as shown. (D) Depiction of the 
effective brick lattice. Dashed rectangles outline a unit cell of the 
lattice. 


energy scales. An analysis of the energy scales of a physically 
realistic system indicates that the Ising topological order could 
have energy gaps on the order of a few percent of the charg¬ 
ing energy of the Josephson junctions of the system; given 
present-day materials and technology, and the constraints on 
the required parameter regimes, we estimate the possibility of 
energy gaps of up to several Kelvin. 


Second, we show how short overpasses between neighbor¬ 
ing chains, which are feasible with current nanofabrication 
technology, can be used to create two effectively independent 
Ising phases, referred to as an Ising x Ising state. Changing 
the connectivity of the network by creating a lattice disloca¬ 
tion allows the creation of a genon; this effectively realizes a 
twist defect that couples the two layers together. Finally, the 
genons can be effectively braided with minimal to no phys¬ 
ical movement of them, by tuning the effective interactions 
between them. 


Realizing the Kitaev Model 

We begin by providing a physical realization of a Kitaev 
model ES, which can be described by the following Hamil¬ 
tonian with spin- 1/2 degrees of freedom on each site r of a 
brick lattice CH: 

r R 

where for a = x, ^ are taken to be the Pauli matrices. 
We take the unit cell of the brick lattice to be two vertically 
separated neighboring spins; ^ two-spin 

unit cell and R refers to the top-most spin within the unit cell. 

i^K is most naturally solved by expressing the spins in terms 
of Majorana operators 7 -^ as = 27 ^ 7 ^, together with a 
gauge constraint Yla=x,y,z,t^^ = ^ ESI- can thus be 
rewritten as 

Hk = J2 JyxlplhArfr+x + E 

r R 

( 2 ) 

where we have made use of the gauge constraint in writing the 
Jz term. The gauge constraint allows us to separate the Majo- 
ranas into a set of Z 2 gauge fields Uxy{R ^ x/2) = 

and Uzz{R + ^/ 2 ) = ^ 7 ^ 7 ^^^’ which commute with to¬ 
gether with Majorana modes 7 ^ that couple to these Z 2 gauge 
fields. The non-Abelian Ising phase corresponds to the regime 
where 7 ^ forms a topological superconductor, such that the 
Z 2 vortices of u localize MZMs. Since the Z 2 vortices are 

deconfined finite-energy excitations in the Ising phase, their 
topological twist is well-defined and can be exploited 
for a topologically protected ir/S phase gate. 

The Majorana solution of i^K suggests that it might be 
physically realizable in a system where each spin is rep¬ 
resented using a pair of proximity-induced superconducting 
nanowires with four MZMs 7 ^’^’^’^, as shown in Fig. [^. 
The semiconducting wires (such as InAs or InSb wires) can 
be either grown 1541 or lithographically defined on 2D sys¬ 
tems na, and the superconductor thin films (such as Al, and 
perhaps other superconductors such as Nb and NbTiN) can 
now be epitaxially grown with exceptional interface qualities 
csiiiiEa. As we show below, such a physical realization 
of i^K is indeed possible; we first provide a proposal that 
physically implements the effective spins using the physics 
of charging energies, and subsequently we provide a second 
proposal which utilizes quantum phase slips. 

In our first proposal (see Fig. [^, an overall charging en¬ 
ergy for the pair of islands, each of which is controlled by 
a capacitance Cg to a gate placed at a voltage Vg, is used to 
effectively generate the gauge constraint by constraining the 
total charge of the pair of islands. For the four MZMs from 
each A, B nanowire pair to be coherent with each other (in a 
sense which will be made more precise below), it is necessary 
to retain some phase coherence between the neighboring is¬ 
lands that comprise a single effective spin. This is obtained 
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by connecting the islands A and B with a Josephson junction, 
with Josephson coupling Ej and capacitance Cj. 

The effective Hamiltonian describing the Majorana and 
phase degrees of freedom for the system in Fig. [JA) can be 
written as 

Hss= E 

3=A,B 

- Ej C0S{(PA - </>b) + ^ E Qi^r/Qj- (3) 

i,j=A,B 

Here is the BdG Hamiltonian for the 

nano wire on the jth island, where |Aoj| is the proximity- 
induced superconducting gap on the jth nanowire (at zero 
magnetic field). Qj = e{—2id^. + Nj — rioffj) is the excess 
charge on the jth superconducting island - nanowire combina¬ 
tion. —id^. represents the number of Cooper pairs on the jth 
superconducting island, Nj = f is the total number of 
electrons on the jth nanowire, and rioffj is the remaining off¬ 
set charge on the jth island, which can be tuned continuously 
with the gate voltage Vgj . The capacitance matrix is given by 
Cij = {Cg E 2Cj)Sij — Cj, for j = A, B. 

In order to decouple the fermions 'ipj in Ffedc from the 
phase fluctuations ipj, we perform a unitary transformation 

u = UHssW = 

Ej ^^BdolAoj, V’j,■*/’j] + H+ + H_, with 

H- = Ec- - Ej cos((/?_), 

H+=Ec+(n'+ + 'A1±ZA^^ , ( 4 ) 

where we have defined if- = ip a — ^b, Ec- = 4e ^/{Cg -h 

2Cj), Ec+ = ICg, rioffi = noffA ± = 

uma i bmb, and = —id^^j2. For wires A^B in the 
topological superconducting phase, at energies below the sin¬ 
gle particle gap Aj, i^Bdc creates essentially decoupled Ma¬ 
jorana zero modes 7 -^ which affect the phase dynamics only 
through the occupation numbers uma — (1 + ^ 7 ^ 7 ^)/ 2 , 
bmb = (1+i 7 ^ 7 ^)/ 2 . In order to allow the MZMs to remain 
free except for a constraint on the total fermion parity, we con¬ 
sider tuning the gate voltages so that noff+ = 2m E 1, where 
m is an integer, and noff- = 0. The ground state of the system 
is then two-fold degenerate, with N'_^ = m/ 2 , nM+ = 1 , and 
um- = ±1. There is an energy gap on the order of Ec+ to 
violating the gauge constraint by changing the total charge of 
the system, and, for Ej < Ec-, a gap of order Ec- to ex¬ 
cited states of that are related to fiuctuations of the relative 
phase if-. 

For energy scales below Ec±, the system can therefore be 
described as an effective spin-1/2 system, with = um- = 
±1. Tuning noff- slightly away from zero acts like a Zeeman 
field for this effective spin degree of freedom, giving an effec¬ 
tive Hamiltonian iifeff,ss = hzS^, with hz oc noff-. This is 
equivalent to tunneling terms between the Majorana modes 7 ^ 


and 7^. Effective Zeeman fields hxS^, hyS^ can also be in¬ 
duced by allowing electrons to tunnel between the MZMs 7^, 
7^ and 7^, 7^ (see Supplemental Materials). We emphasize 
that the MZMs 7^ are not exactly equivalent to the 7^ used 
in Eq. because the two effective spin states in this setup 
differ not only in the Majorana occupation numbers um-, but 
also in the wave function of ip - (see Supplemental Material). 
We can think of 7 as corresponding to 7, dressed with the ip- 
degrees of freedom. 

The next step is to generate the quartic Majorana couplings 
in Eq.j^by coupling the different A — B island pairs together. 
For example, the Jz term in Eq. [^essentially represents a cou¬ 
pling between the occupation numbers nM,A,f and um.Bc+z 
of neighboring SC islands in the lattice and can be realized us¬ 
ing a capacitor Cz- This is shown in Fig. [JC), where we have 
also introduced the labelling 1,..., 4 for the vertically coupled 
islands. For numerical optimization of energy scales, it is use¬ 
ful to also consider a capacitance C'^ between islands 1 and 4, 
which is not shown explicitly in Fig. 

To estimate the resulting Jz coupling, we consider a de¬ 
tailed model for two vertically coupled effective spins, con¬ 
sisting of four vertically separated islands (see Fig. m 
which is described by a Hamiltonian 

H 2 s =Hi 2 + H 34 + ^1234- (5) 

Hi 2 and 7^34 are the Hamiltonians for the iso¬ 
lated units of the form of Eq. while H 1234 = 

E<ti,<t 2 =± * 512,^1 Q 34 ,<72^ai<T 2 capacitively couples the 
two effective spins. As expected H 1234 couples the differ¬ 
ences Qij^± = Qi E Qj of the charges Qj on the islands. 
Acricrs = i {Ci^ E cr 2 Ci^ E criC^s + ^ 2 ^) are related to the 
four-island capacitance matrix (see Supplemental Materials). 

The term H 1234 in Eq. generates a coupling Jz between 
the effective spins 1 — 2 and 3 — 4 by coupling the charges 
Qj on the various islands. For small coupling capacitances 
Cz, C'z, this can be estimated perturbatively. The limits of 
validity of the perturbative estimate can be checked by a di¬ 
rect numerical calculation of the spectrum of H 2 S, which 
we have performed and presented in the Supplemental Ma¬ 
terial. An example of a suitable parameter regime is for 
Cz = C'z = O.bCg, Cj = l.bCg, and Ej = OAbe^/Cj. 
In this case, we find Jz ~ 0.02e^/Cj, while the gap to all 
other states in the system is F^gauge ~ lOJ^. Thus the gauge 
constraint is implemented effectively through a large energy 
penalty F^gauge. and the system is well-described at low en¬ 
ergies by the effective spin model (or, equivalently, the con¬ 
strained Majorana model). 

The Jxy terms in Eq. [^involve coupling MZMs in the hor¬ 
izontal direction. This quartic Majorana coupling can be ob¬ 
tained from single electron tunneling processes between the 
MZMs through (normal) semiconductor wires that run hori¬ 
zontally, as shown in Eig. [Tp. The electron tunneling ampli¬ 
tudes tA and Ib can also be controlled with a gate voltage. 
The resulting Hamiltonian for the full 2D system shown in 
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Fig. [Tp,D is then 


A 


B 


H2D — ^ ^2s,R ^ 

R 

H\xxn — ^ “1“ ^-B'0y^^'0cc,r+:E H“ H.C.\^ (6) 

f 

where ^ is the Hamiltonian for the two-spin unit cell at 

R, given by ^ above. The single electron tunnelings Ia, 
Ib violate the gauge constraint, which is related to fermion 
parity of the single effective spin, and thus induce an energy 
penalty on the order of F^gauge- We consider the limit where 
IaRb ^ ^gauge, SO that Fftun Can be treated perturbatively 
around the decoupled unit cell limit. Assuming further that 
Ia "C A^, Ib ^ A^, where Aj is the single-particle gap on 
the jth superconducting nano wire, we can replace after 
the unitary transformation U, by the MZMs: 

UHa,rU = (7) 

where we have set after the unitary trans¬ 

formation U, the c-number is the wave function of 
the MZM, and FpA,p = FpB,r = hpll are the 

fermion parities of the A and B islands of site r. It is 
useful to define is = tBUy^^Ux^r+x, 

so that after the unitary transformation by U, Fftun = 

H.c.]. Treating Fftun perturbatively around the decoupled unit 
cell limit, we obtain an effective Hamiltonian 


H,ff = 


tAtB 

^gauge 




Cab are constants that depend on parameters of the model, 
in particular the angle 0 = ArgitABB)' « 0 , 7 r and 

Ej « Ecj = find Cyx > Cxx,Cxy,Cyy. There¬ 

fore, up to n^ligible corrections, we arrive at the effective 


Hamiltonian ( l| above, with Jyx = 
the angle 0 can 


Cyx. Physically, 

-t^gauge ^ 

~ 3 e tuned by the applied magnetic fiux piercing 
the loop defined by the tunneling paths IaRb and also the an¬ 
gle between the Zeeman field and the Rashba splitting of the 
wire. 

In the Supplemental Materials we provide detailed numer¬ 
ical estimates for the parameters of the effective spin model. 
For the parameter choice Cz = C'^ = 0.5C^, Cj = 1.5C^, 


Ej = 0.45e^/Cj described above, and jt^l, \ti 


: 0.2K 


0 ^ O^TT, we find that \cy 


1.75, |Ca;; 7 ;| — 0.3, Cx 


gauge 

0 , 


and thus Jyx ~ 0.016e^/Cj. Combined with the above esti¬ 
mate Jz ~ 0 .02e‘^ jCj, we see that the energy scales Jyx, Jz 
are on the order of a few percent of the charging energy of 
the Josephson junctions. For Al-InAs-Al Josephson junctions, 
one can achieve (561 ~ 3 — 4 K ^ Ecjl2, implying that 

Jyx 5 Jz would have energy scales on the order of several hun¬ 
dred milli-Kelvin. Nb-InAs-Nb Josephson junctions would in 



FIG. 2: Decorated brick lattice, to realize variant of the Ki- 

taev model (m (A) A single plaquette shown made of the 
superconductor-semiconductor system. (B) Effective lattice model. 
Each unit cell now consists of 6 effective sites of the lattice. 


principle be able to support energy scales on the order of 7-8 
times larger, with Jz,xy on the order of several Kelvin, due the 
correspondingly larger superconducting gap of Nb. 

The phase diagram of Eq. contains two phases (531 : 

a phase where the fermions 7 t form a trivial insulating phase 
when Jz ^ Jyx, and a phase where the fermions 7 ^ are gap¬ 
less with a Dirac-like node when Jz ~ Jyx- Both phases 
have an Abelian topological sector associated with the Z 2 
gauge fields Uxy^z- ft is possible to open a topological gap 
for the Dirac node, and thus realize the non-Abelian Ising 
phase, by breaking the effective time-reversal symmetry of Q 
with a Zeeman field ^f^y=x,y,z the implementation 

of which was described above. 

A potentially more optimal approach to inducing the non- 
Abelian Ising phase is to use a modified structure where each 
point of the brick lattice of Fig. is expanded into three 
points, with the couplings as shown in Fig. [^ (57l . The ground 
state on this lattice spontaneously breaks the effective time- 
reversal symmetry of Q and gaps out the Dirac nodes in the 
regime where Jz ^ Jxy to open a topological gap on the order 
of Jz ^ Jxy 

In our system, the small perturbations hz, Cxy, Cxx, Cyy can 
be used to controllably tune the sign of the spontaneous time- 
reversal symmetry breaking and thus control whether the sys¬ 
tem enters the Ising phase or its time-reversed conjugate, de¬ 
noted Ising. 


Ising X Ising phase and genons 

A crucial feature of the physical setup that we have pro¬ 
posed is that the vertical couplings between neighboring spins 
only involve capacitances (or Josephson junctions as de¬ 
scribed in the quantum phase slip based implementation be¬ 
low). This means that once a single copy of the model is re¬ 
alized, it is straightforward to realize two effectively indepen¬ 
dent copies of the model by creating short overpasses. Specif¬ 
ically, this can be done as shown in Fig. |^by fabricating the 
superconducting wires that run in the vertical direction to pass 
over one pair of nanowires and to couple capacitively to the 
next chain over in the vertical direction. 

The Ising phase contains three topologically distinct classes 
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FIG. 3: Two effectively decoupled copies of the Kitaev model can 
be created with short overpasses where the vertical superconduct¬ 
ing wires skip over one chain and couple to the next chain. Genons 
(green circles) can be created at the endpoint of a branch cut along 
which the vertical couplings connect the two copies of the model 
together. 


A B 



FIG. 4: Genons can be braided without moving them by sequentially 
projecting different pairs of genons onto specific fusion channels, as 
depicted by the gray ellipses in A - D. This can be implemented by 
tuning the interactions between the genons by decreasing the quasi¬ 
particle gap along various paths. 


of quasiparticle excitations, labelled as I, t/), and a. The Ising 
X Ising phase contains nine topologically distinct classes of 
quasiparticles, which we label as (a, 6), for a, 6 = I, t/;, cr. 

A genon, which we label Xi in the Ising x Ising phase is 
a defect in the capacitive couplings between vertically sepa¬ 
rated chains, associated with the endpoint of a branch cut that 
effectively glues the two copies to each other (see Fig. 
ca. This defect in the lattice configuration of the supercon¬ 
ducting islands is not a quasiparticle excitation of the system, 
but rather an extrinsically imposed defect with projective non- 
Abelian statistics |[3T]|58l 1761 . Xi has quantum dimension 2, 
and possesses the following fusion rules 1771 : 

XiY. Xi = (I, I) + + (cr, cr). ( 9 ) 


the system to effectively project the genons 5 and 6 onto the 
fusion channel 656, then the genons 5 and 3 onto the fusion 
channel 635, the genons 5 and 2 onto the fusion channel 625, 
and finally again the genons 5 and 6 onto the fusion channel 
650 (see Fig. 22 ). We will asssume here for simplicity that 
^56 = ^56- If the genons 5 and 6 are created out of the vac¬ 
uum, then it will in fact be natural to have 656 = 655 = (I, I). 
As long as 650, 635, 625 are Abelian, i.e. equal to either (I, I) 
or (V^, V^), then the results of OTl [ 5 ^ imply that the matrix 
obtained for a double braid (i.e. a full 27 r exchange), is given 
by 


/10 0 \ 

(i?23)' = e'G 0 1 0 , (10) 

\0 0 


In ED, it was shown that the braiding of genons maps to 
Dehn twists of the Ising state on a high genus surface, which is 
known (381 [40l to provide a topologically protected tt/S phase 
gate. The protocol for implementing the tt/S phase gate us¬ 
ing genon braiding was described in ED In the present sys¬ 
tem, the braiding of genons is complicated by the fact that 
it is difficult to continuously modify the physical location of 
the genons to execute a braid loop in real space. Fortunately, 
this is not necessary, as the braiding of the genons can be 
implemented through a different interaction-based approach, 
without moving the genons, as described for general any on 
systems in (Wll . To do this, we require that it be possible to 
project the joint fusion channel of any pair of genons into ei¬ 
ther the (I, I) channel or the channel. This can be done 

by adiabatically tuning the effective interactions between the 
genons, similar to proposals for braiding MZMs in nanowire 
networks (601 [Ml . 

In order to implement the tt/S phase gate, we wish to start 
with two pairs of genons, labelled 1 , 4 , and have the abil¬ 
ity to braid genons 2 and 3 . To do this, we require an ancil¬ 
lary pair of genons, labelled 5 and 6. The braiding process is 
then established by adiabatically changing the Hamiltonian of 


where is an undetermined, non-topological phase. In other 
words, the state obtains a relative phase of if the fusion 
channel is (cr, cr) as compared with (I, I), or (7/^, 7/^). 

When two genons are separated by a finite distance L, the 
effective Hamiltonian in the degenerate subspace spanned by 
the genons obtains non-local Wilson loop operators: 

-f^genon = “1“ ^cr^^(a,a) “1“ H.C. (H) 

W(^a,a) describes the exchange of a (a, a) particle between 
the two genons, which equivalently corresponds to a (I, a) or 
(a, I) particle encircling the pair of genons. ta oc 
for a = are the tunneling amplitudes, with being 

the energy gap for the a quasiparticles, and Va their veloc¬ 
ity. When a ( 1 , a) quasiparticle encircles a topological charge 
( 6 , 6 ), it acquires a phase Sah/Sih, where S is the modular S 
matrix of the Ising phase El- As we show in the Supplemental 
Materials, this implies that for our purposes, we only need to 
ensure that ^ 0 and, in the case where \ta\ < |, we must 

have < 0 . 

The tunneling amplitudes ter can be tuned physically by 
tuning the parameters of the model, such as the electron tun¬ 
neling amplitudes Ia, ts, the capacitances Cz, and the gate 
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voltages Vgj. Therefore, to tune the interactions between two 
desired genons, we tune the parameters of the model in or¬ 
der to decrease the energy gap to the quasiparticle excitations 
along the path that connects them. A more detailed study of 
this will be left for future work. 

We note that it is also possible to implement effectively the 
same physics by using instead the Ising x Ising state. In this 
case, the genons are replaced by holes with gapped bound¬ 
aries, and the topological charge projections are implemented 
along various open lines that connect the different gapped 
boundaries. A detailed discussion of this variation is pre¬ 
sented in the Supplemental Materials. 
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Quantum Phase Slip Limit 


An alternative architecture is also possible, if we replace the 
purely capacitive coupling Cz by a Josephson junction, with 
Josephson coupling Ej^ and capacitance Cz • In this case we 
consider the limit where the Josephson energies Ej, Ej^ are 
much larger than the charging energies leading to a 

state with long range phase coherence. In the limit where the 
charging energies are ignored, the system has a large degen¬ 
eracy due to the MZMs. A small charging energy induces 
quantum phase slips of the superconducting phase of the is¬ 
lands; the amplitude of the quantum phase slips depends on 
the occupation of the MZMs on the superconducting islands, 
thus inducing an effective Hamiltonian in the space of states 
spanned by the MZMs. The effective Hamiltonian takes the 
form: Ff 2 D = + i ^2 + ^tun, where Hi consists of single¬ 

island phase slips: 


Hi = 


( 12 ) 


and H 2 consists of double island phase slips: 




-5ICz“7g757B_.7' 


AB^^x y t 

R ^R-z ^R-z 


(13) 


Fftun is the same as in Eqn. and describes electron tunnel¬ 
ing in the horizontal direction. The single island phase slips 
are modulated by the offset charge: ("i ex cos(7rnoff,r,j) and 
can therefore be tuned to zero using the gate voltages. Dou¬ 
ble island phase slips that are not included in H 2 can be ig¬ 
nored in this limit, as can phase slips that involve more than 
two islands, as they are exponentially suppressed. We wish 
to choose parameters to operate in the limit Cz^ ^ 

0^. In this case, effectively imposes the gauge con¬ 
straint 7 ^ 71757 ^ = 1 for states with energies much less 
than 0^. Each site can therefore be described by a spin- 
1/2 degree of freedom. The term involving 0^ acts like a 
coupling As before, in this limit iiftun can be 

treated perturbatively, and gives rise to the desired coupling 

I qy qx 
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FIG. 5: (a) A pair of A, B superconducting islands, coupled together with a Josephson junction, denoted by the x. The red dots represent 
Majorana zero modes, labelled 7 ^, 7 ^, 7 ^, 7 ^. The Majorana zero modes are localized at the end of the semiconducting nanowires, which are 
shown in black. Each A and B island has a capacitance Cg to a gate voltage Vg, which is then grounded, (b) An effective circuit for the pair of 
A, B islands, indicating the capacitances and Josephson junctions, (c) Grey wires indicate normal (non-superconducting) semiconductor wire. 
Tuning the voltage on the semiconductor wires allows tuning the electron tunneling amplitudes tx and ty between the Majorana zero modes. 


Supplemental Material 
Charging Energy Based Implementation 

We begin by describing an implementation where charging energies are the dominant energies in the system. We proceed 
by incrementally increasing the complexity of our analysis, by first describing the physics of a single effective spin, then two 
vertically coupled spins that will form the unit cell of the Kitaev model, and subsequently the horizontal couplings that will link 
all of the two-spin unit cells into the full effective spin model. 


Single spin 

Let us consider the configuration shown in Fig. |^a), which consists of two superconducting islands, labelled A and B, each 
of which is proximity coupled to a Majorana nanowire. Each superconducting island is separated by a capacitance Cg to a gate 
voltage Vg. The A and B islands are coupled together through a Josephson junction, with Josephson coupling Ej and junction 
capacitance Cj. The effective Hamiltonian for this system is 

Hss= E + \ E (14) 

j=A,B i,j=A,B 

Here, ipj for j = A^B is the superconducting phase on the A and B islands, is the BdG Hamiltonian for 

the nano wire on the jth island, where | Aoj | is the proximity-induced superconducting gap on the jth nano wire at zero magnetic 
field. Qj is the excess charge on the jth superconducting island - nanowire combination; it can be written as: 


Qj = e{-2id^. + Nj - nofy), (15) 

where represents the number of Cooper pairs on the jth superconducting island, Nj = f ip^jipj is the total number of 

electrons on the jth nanowire, and riogj is the remaining offset charge on the jth island, which can be tuned continuously with 
the gate voltage Vg. 

The capacitance matrix is given by 


fCg + Cj -Cj \ 
\ -Cj Cg+CjJ 


( 16 ) 
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The charging energy term can be rewritten in terms of the total and relative charges on the A and B islands: 

2 ^ ^ 4 Cg + 2Cj ~ 

2 2 

= ^(-idvA - idvB + (^+ - 'rios+)/2f + ^ ~ (17) 

where we have defined 


N± = Na ± Nb, 

^off± = ^offA ± ( 18 ) 


The BdG Hamiltonian for the nano wire is given by 


iJEdolAe*"^, J dx[A{x) + iaayd^ + qubB • ^p{x) 

+ (Ae*‘^V’tV’| + 77-c.)], 


(19) 


where 


. Here we have taken x to be the coordinate along the wire and L is the length of the wire, a is the Rasha 


spin-orbit coupling, fi is the chemical potential, and m* is the effective mass of the electrons in the nanowire, B is the magnetic 
field and gjj.B\B\ is the Zeeman energy. 


It is now useful to perform a unitary transformation U = e 




in order to decouple the phase (^j from 


the fermions pjj in i^Bdc- Here, umj = 0,1 is the occupation number of the pair of Majorana zero modes on wire j. It is given 
in terms of the Majorana zero modes as 


riMA = (1 +n^7*)/2, 

nMB = (l+n"7^)/2. (20) 


Under this transformation, the charge Qj transforms as: 

Qj — U^QjU — 6( “h ^Mj ^offj )• 

Thus, taking Hgs U^HssU, we obtain 


=U^H,,U = Y,HB6G[Aoj,Aj,i’j] + H+ + H_, 


4e^ 

H+ — + 


C. 


H- = 


g 

4e^ 


Bm+ — Boff_+ 


Cn + 2Cj 


-idu: 


um- — rioff- 


- Ej cos(/p-). 


Here, we have defined the combinations: 


nM± = nMA ± nMB, 

‘P± =‘PA=t: ^pB, 

With this definition, 

[id^^{x),ifip,{x')\ = i5pp>S{x - x'), p,p' = ± 

It is also useful to define 


( 21 ) 


( 22 ) 


(23) 

(24) 


Nj = -idpj, j = A,B 
N:^ = {N'^+N's)/2 

Q'± = Q'a a Q'b- (25) 
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With these definitions, we see that the compactification of Lp- is: 

((/?+,(/?_) ^ ((/?+ + 27r, if - + 27r) ^ ((/?+ + 27r, (p_ - 27r). (26) 

While (/?+ and (/?_ are formally decoupled in the Hamiltonian, they are coupled through their boundary conditions. The ground 
state of H'gg can now be written as 

1 

- / dip+dip_^riMA,nMB{‘f+^‘f-)\‘f+,‘f-^'n'MA,nMB) (2.1) 

^ J- 27 T 


Here, {uma^ i^mb) is the state of the Majorana zero modes, and (/?_) is the state for the phase degrees of freedom. Impor¬ 
tantly, the wave function of (/?+, (f- itself does depend on the values of uma, ^mb- Since the Hamiltonians for and H- are 
decoupled, we can immediately write the ground state wave function for ip_: 

'>PnMA,nMB{<P+,V’-) = ( 28 ) 


where fuM- (^-) i^ ground state wave function for which is peaked at if- = 0. Importantly, because of the com¬ 
pactification conditions on (p±, we see that if N'_^ = + ^b)/‘^ i^ integer (half-integer), then fuM- (^-) must be periodic 

(antiperiodic) in (/?_ (p_ ^2 tt. 

For energies much less than the single-particle gap A on the nanowire, we can ignore the excited single particle states associ¬ 
ated with FfBdG, and we can describe the system by the following effective Hamiltonian: 

^M+ 


Let us define 




— -C- 


4e^ 


2 Cj 


-id^- -h 


Um- — ^off,- 


Ej cos((/?_) 


(29) 


46^ 

Ec+ = 7^, 


Ec- 


4e^ 


Cq -f 2 Cj 


(30) 


We see that the effect of is to fix the total charge N'_^ on the pair of A, B islands to a fixed value, and gives an energy cost 
of Ec+ to increase the charge by one unit. If we set 


noff+ = 2m + 1, m e Z, (31) 

then the ground state of the system will be given by 

N'_^ = m/2, nM+ = 1- (32) 

Therefore the system will have two lowest energy states, associated with um- = ±1. Thus we can define an effective spin 
degree of freedom: 


= riM- = El. 


(33) 


We will denote these two states as \S^)\ 


I gz I _ gz 

IS'^) oc / d(p+dip-fsz{ip_)\ip+,ip_,nMA = —^—,nMB = —;;—), 


(34) 


where we have chosen N'_^ = m/2 = 0 for simplicity. Note that \S^ = ±1) differ both in the value of um- = E and also the 
wave function fuM- (^-)- future reference, it will also be useful to define the state 

_ I gz I _ gz 

oc J dip+dip-f-s^iip-)\^+iV-^'nMA = —^— ,nMB = —^—), (35) 

which has the opposite wave function /_£;* {(/?_) as compared with 15^). 

The effective Hamiltonian in the two-dimensional space = ±1) is given (up to an overall constant) by 


H,n = h,S^. 


( 36 ) 


The value of depends on parameters in iJeff, as described below. 
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EI\Ec Ej = 0 



FIG. 6: Plot of energy spectrum of H-, in units of AEc. We have set Ej — 0. Green and red curves are the ground state and the first excited 
state energies, respectively, for the case — tim- — 1, as a function of the offset charge rioff-. Blue and yellow curves are the ground state 
and first excited energies, respectively, for the case — um- = — 1. We can see that at noff- = 0, there is a degeneracy between = ±1, 
with a gap of order Ec to all other states. Non-zero rioff- acts like a Zeeman field that splits the energies of the two spin states. 


E/AEc Ej/AEc = 0.2 



FIG. 7: Same plot as in Fig. but for Ej — ^.2(AEc). 


Numerical Solution 


The Hamiltonian (see Eqn. (22)). The three terms, commute with each other and 

can be separately solved. As discussed above, i^Bdc is gapped for energies below A, aside from the zero energy states arising 
from the Majorana zero modes, while has a gap of F^c'+- H is useful to solve numerically, for the different Majorana 
occupation numbers. 


In Fig. I^and Fig. we plot the energy spectra for the four lowest energy states of as a function of the offset charge 
noff,-, for the two different values of the Majorana occupation numbers = um- = ±1- To connect with some standard 
notation in the literature for the well-known Hamiltonian H -, it will be useful to define 


Ec = Ec-!^ 


CgE2Cj' 


(37) 


We see that for energies much smaller than Ec, the system simply consists of the two states \S^ = ±1). These are degenerate 
when noff- = 0, and acquire a small splitting when noff- ^ 0. 
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Analytical Solution 

The Hamiltonian can also be fully solved analytically through the use of Mathieu functions (see, e.g., Ref. |63l). Here we 
will provide this solution for reference. We find that 

hz = Ec{ay^{—Ejl2Ec) - a^_{-Ejl2Ec))l2 (38) 

where a^{(i) is Mathieu’s characteristic value, and 

= 2[-ng± + /c(0,n^±)], (39) 


ng± = (±1 -noff-)/4 + 


m mod 2 


/c(0,n^±) = int(n^±) ^ (int(2n^± +//2) mod 2). 


i=±i 


(40) 


int(x) rounds x to the nearest integer. The average charge Q'^ on the B island is given by: 

{S^\Q'b\S^) = = -e{S^\ - 2id^_ + (S^ - no«,-)/2\S^) 

e {s^idH'jsy e dh. 


Ec driofi- 


Ec driQii- 


Evaluating the partial derivative: 


- a',_{-Ejl2Ec)), 


(41) 


(42) 


where 


, da^{x) 


(43) 


Note we have assumed that dkjdng^ = 0, which is true except for certain fine-tuned values of 


Effective and S'^ terms 

Above we showed that tuning noff- away from zero effectively acts like a Zeeman field in the direction. Zeeman fields in 
the and direction can also be generated, by allowing electron tunneling, with amplitude tx and ty, through the semicon¬ 
ducting wires as shown in Fig. [^. Consider the following electron tunneling perturbations to Hss- 

SH = + H.C. (44) 

After the unitary transformation U, SH changes: 

5H' = U^SHU = + H.C., (45) 


where 

= u'^'ipyU = 

t/,' = (46) 

where 

FpA = 

FpB = 


( 47 ) 







13 


are the fermion parities of the A and B islands, respectively. Assuming the regime 

txAy ^ (48) 

where A is the single-particle gap in the semiconducting nanowire, we can write the electron operators at low energies in terms 
of the Majorana zero modes: 

(49) 

where a = z, t, and Ua are complex numbers (whose magnitude is order unity) that depend on microscopic details. Thus, 
we obtain: 

Recall that nM+ = 1, and that uma = (1 + FpA)l‘2, umb = (1 + FpB)l‘2, which implies htat 

FpA + FpB = 0. (51) 

It is useful to define 

tx ~ 

iy = tyulut (52) 

Thus, we get 

= l^^-i^B{l-FpA)/2+iy:>A{l-FpA)/2^x^t :^^^-iy:>B{l-FpA)/2+iy:>A{l-FpA)/2^y 

= ^ jj ^ ^ ( 53 ) 

where we have also commuted 7 ^, 7 ^ through the exponential term. Note further that 5'^ = um- = ^ma — ^mb = {BpA — 
Fpb)I = FpA- The above can then be rewritten as 

SH' = (47"7* + + H.c., 

= j + iyAl') + H.C., 

=i7"7*Im ( 4(1 + e*^-) - 4(1 - 4 ^-)) 

+ ( 4(1 + e"'^- ) - 4(1 - )) (54) 

Thus, we have 

SH' = (55) 

with 

= Im (4(1 + - 4(1 “ e*‘^“)) 

Gy = Im (4(1 + - 4(1 - e*‘^“)) (56) 

For txAy ^ we can treat 6H' perturbatively around Hss- Thus, we get an effective Hamiltonian Ffgff, such that 

{m\H,ff\n) = {m\Hss^SH'\n) 

= SmnEm + (57) 

where \m) are the normalized eigenstates of the unperturbed Hamiltonian Hgs- Then we can write the effective Hamiltonian in 
the low energy spin space as 

Heff = '^Em\'ni){m\ + '^{m\6H'\n)\m){n\ 

m m,n 

= h.S"' + + hyS\ 


( 58 ) 
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FIG. 8: Two vertically coupled spins, consisting of 4 superconducting islands. Each island is labelled 1,...,4 as shown. Cz is a capacitor 
connecting islands 2 and 3. In order to optimize energy scales, we will also consider a capacitance C' connecting islands 1 and 4, which we 
have not explicitly shown. 


where 


h, = Es^=i - Es^=-i + {S^ = 1\SH'\S^ = 1) - {S^ = -1\SH'\S^ = -1), 
ft, =Re[(5^ = -l|5if'|5^ = l)], 
hy=lm[{S^ = -1\SH'\S^ = 1)]. 


Recall that the two spin states of interest, jS'^) are defined as in Eq. (34). Thus: 


= 0 , 

{S^\6H'\ - S^) = (5^|a,z7"7* + - 5^) 

= {S^\a,\W)-iS^S^\ay\^), 


where \S^) is defined in Eq. (35). Therefore, we find: 


/i, =Re[(-l|a,|-l)+i(-l|a^|-l)] 

hy = Im[(-l|a^| - 1) + i{-l\ay\ - 1)] 


(59) 


(60) 


(61) 


Two spin unit cell 

Let us now consider a pair of vertically separated effective spins, which will form the unit cell for our brick lattice Kitaev 
model. This consists of two pairs of A and B islands, as shown in Fig. We wish to consider a capacitive coupling Cz, as 
shown in Fig. In this analysis we will label the islands 1,4 as shown. We will also consider a capacitance C'^, purely for 
subsequent numerical optimization of energy scales, between islands 1 and 4, though this is not explicitly shown in Fig. The 
Hamiltonian for such a two-spin system is given by 

^ ^ 1 . 

H 2 s = + 2 51- Ej{cos{(fii - (^ 2 ) + cos((y?3 - (^ 4 )). (62) 

i=l ij 

The charges Qj (after the unitary tranformation discussed in the previous subsection) are 



Qj — 2e 


2 


(63) 
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Note we omit the primed superscripts in the preceding equation and throughout the rest of the discussion. The capacitance 
matrix C is now a 4 x 4 matrix: 


(Cg + Cj + C'z -Cj 0 

-Cj Cg + Cj + Cz -Cz 
0 -Cz Cg + Cj + Cz 

-C'z 0 -Cj 


C'z ^ 
0 

-Cj 

Cg+Cj+C'zj 


It is useful to write H 2 s as 


H2s — Hi2 + Hs4 + -ffl234, 


(64) 


(65) 


where 

2 ^ 2 ^ 

Hi 2 = + 2 51 ~ -E’jCOs((y?i - ip 2 ) 

j = l i,j=l 

2 

= ^ -ffBdG[Aoj, t/’j, tAj] - Ej C0S((^12,-) 

i=i 

+ q? 2,+J(6?riV2 + C 22 V 2 + cr2') + q? 2,-J(6?n/2 + C 22 V 2 - cr^i) + gi2,+Qi2,-- c^,^) 

4 ^ 4 

^34 = ^ ^BdG[^0i: + 2 S ~ “^ 4 ) 

j =3 i,j =3 

4 

= ^i^BdG[Aoj,' 0 ],' 0 j] - Ej COs((p 34 -) 

J =3 

+ Qi4,+ l(6?33V2 + C 44 V 2 + C 3 -/) + Qi4,-l(6?33V2 + C 44 V 2 - C 3 -/) + g34,+Q34,-l(C3-3l - C 4 -/) 

2 4 

^1234 = 

i=l j =3 

= E *312 ,ai Q 34 ,(J 2 ^<J 1<J2 (66) 

cri,cr2=± 


where we have defined Qij,± = Qi E Qj and 


^(Ti(T 2 — ^(^13 + H- CF 1 C 23 + C 24 ) (67) 

The terms H 12 and 7^34 are just the Hamiltonians for a single effective spin, which was analyzed in the previous section. We 
label these spins by and S^_-. r and f—z label the two different effective sites, as shown in Fig. i7i234, then, couples the 
two effective spins. 

Analytical treatment 

We now wish to treat H 1234 perturbatively around the decoupled limit i 7 i 2 + 7734 - This is valid if 

< Ec. (68) 

To lowest order in perturbation theory, we can replace H 1234 with the effective Hamiltonian Ei^ay. 

(rSSrSS_,|77eff,,|rSSrSS_,) = | 77l234 | ) 


( 69 ) 
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Cj = l.SCg, Cz = 0,Ej = 0 E/Ecj 




1 1 E 


2.0 

Cz/Cg 


Cj = l.SCgXz = 0,Ej = 0 


Cz/Cg 


FIG. 9: Left panel: Plot of energy spectrum of H 2 s (ignoring the excited state spectrum of for the lowest energy states, as a function 

of the capacitance Cz- The other parameters are set to Cj = C'z = 0, = 0. The lowest energy curve is doubly degenerate, and 

is associated with the states S^_^) = (1,1),(—1,—1). The next excited state, whose energy difference with the lowest energy curve 
defines Jz, is also doubly degenerate and associated with the spin states = (1,-1), (—1,1). The next excited state lies outside 

of the effective ’’spin” subspace that we are interested in, and we define the gap to these excited states as F^gauge- The notation F^gauge is used 
because states with energies E > F^gauge can violate the ’’gauge” constraint = 1 that is required for the Kitaev spin model. Right 

panel: Plot of Jz and F^gauge- 


E/Ec, 


Cj = l.SCg, C'z = 0 . 57 , £} = 0 E/Ecj 


Cj - \.SCg,Cz - O.SCgjEj — 0 
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3 * . . . 
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. 
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1.0 
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FIG. 10: Same as in Fig. with different parameters as indicated. 


where we have defined (see eq. 41) 


Jz = A_ 


dhz 


Ec driofi- 


(70) 


The second equality in eq. (69) follows because ((3i2,-) = (Q 34 ,-) = 0, so only the (Qi 2 ,-Q 34 ,-) term remains non-zero. We 
have assumed for simplicity that the two A,B island pairs have the same parameters Ej, Ec, rioff±. 

Therefore, to first order in perturbation theory, the effect of the vertical capacitances Cz, C'^, which couple the two spins at r 
and r — z, is to induce an coupling. 


Numerical Solution 

We can also more comprehensively analyze the two-spin model by employing a numerical solution. In Fig. [^-[TT] we present 
results of such a numerical solution for certain choices of parameters. 
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E/Ecj 


Cj = l.SCg,C^ = O.SCg,Ej = 03Ecj 


0.20 

0.15 

o.ioj; 

0.05 


(S?,S|_^) = (1,-1), (-1,1) 


2.0 

Cz/Cn 


Cj = 1.5Cg,C^ = 0.5Cg,Ej = 0.3Ec 


E/Ec, 




\ 


• • 


Cz/Cg 


FIG. 11: Same as in Fig. with different parameters as indicated. 
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Cj = 13Cg,Cz = ^3Cg,Ej = 0.45Fc 
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1.0 


1.5 


Cz/Cg 


FIG. 12: Same as in Fig. with different parameters as indicated. We see that a particularly optimal point occurs when Cz — 0.5C^. 

Horizontally coupled unit cells: four spins 


Let us now consider horizontally coupling two unit cells, as shown in Fig. 13 We connect two horizontally separated unit 
cells with semiconductor wire, as shown, which allows electrons to tunnel between the end points of the wires, with tunneling 
amplitudes Ia and as shown. The effective Hamiltonian for this system is now 


Has — ^ ^ El2s,I -f^tun) 


(71) 


where H 2 sj is the Hamiltonian for the Jth unit cell, which is given by H 2 s above. i7tun contains the horizontal couplings, as we 
explain below. 

We wish to show that in a suitable parameter regime, at low energies the effective Hamiltonian can be described by the 
following spin model: 

(72) 

with corrections to this effective Hamiltonian being much smaller in energy scale than Jz , Jyx • 

We consider the electron tunneling terms Ia and as shown in Fig[^ This gives rise to an effective tunneling Hamiltonian 
(written in the basis before the unitary transformation U)\ 

77.C.] + “f 77.C.] (73) 

Here, r labels the effective spins, each of which consists of an A and a B island. After the unitary transformation by U, we have 

+ H.C.] 

+ [tB{'>Py,A^'x,r+x + H-C-], 


( 74 ) 
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r+x+z 


r+x 


FIG. 13: Two horizontally coupled unit cells, consisting of 4 spins total. Dashed lines encircle each unit cell. The locations of the spins are r, 
f — z, f-\- X, f-\- X z. tA and ts indicate electron tunneling, as shown. 


where 

(V’a.r)^ = (75) 

where j = A^B depending on whether a = z, t or x, y. Therefore, the tunneling Hamiltonian is, after the unitary transformation: 

Kn + H.C.] 

+ + H.C.], (76) 

We consider the limit where 

IaAb ^ (77) 

where A is the single-particle gap in the nanowire. In this limit, the electron operator can be replaced by 

’ (78) 

where are complex numbers that depend sensitively on microscopic details. Let us also define 

tA = tA{ut)*U^f^+S;, 

(79) 

We therefore write 77/^^ 

T^tun “ 5 

_ J^^^tg-i(l-Fp,A,r)V7Ar/2 + i(l-Fp,A,r+£)¥^Ar+£/2^t^^ 

+ + H.C.] 

— —7(l + ^p,A,'r)¥^Ar/2 + i(l-Fp,A,'r+x)¥^Ar+£/2^t^t^^ 

+ + H.C.] 


(80) 
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These single electron tunneling processes violate the charging energy constraint and are therefore suppressed in the limit 

^ -^gauge: (^1) 

where T^gauge is the energy cost to adding a single electron to the two-spin unit cell. Perturbing in fj/T^gauge, we obtain an 
effective Hamiltonian: 






E, 


gauge 


Expanding, we obtain, up to a constant term. 


^t,eff = — - 




+ ht-2 + ^t;3) 


gauge 

ht-1 = + H.c. 

hi.2 =.‘21 jsJ, ^ ^ ^ ^ 

iWf+xlllf+x + H.c. 

^hhi^hhi + H-c. 

^ —2l^l'^Q~^i^ + FpA,r)^A/2 + i{l-FpA,E+x)^A,r+x/‘^^'ii^-FpB,r)^B,r/‘^-'ii^ + FpB,r+x)^B,r+x/‘^ 

iWr+xlhi'^l + H-C- 
iWr+xlhxll + H.C 

Note that in the limit within which we are working, 

2|?/i?b| 2|f^|^ 2|tB|^ E, 

TP 'TP 'TP ^ -^gauge 

-Hypaupp. -HL/naupp J-jp: 


(82) 

(83) 

(84) 


(85) 


( 86 ) 


(87) 


^gauge -^gauge -^gauge 

Thus, due to the charging energy on each site r, (p^ is highly fluctuating independently on each site. Treating ( [8^ perturbatively 
around the decoupled limit H 2 s, we see that we can set 


Ep-\-^f — Epj^^f -|- FpB^r — 0 ? 

Fp_^f = 2FpA,f = -‘2FpB,r = 25 '^. 

Moreover, we can replace by their expectation values in the ground state manifold of i^ 2 s- 

(^|^t;i,eff|^) = {m\htA\n) = 0 

(^|^t; 2 ,eff|^) = {m\hfp\n) 

ihhxlhhx + H.c.\n)= 0. 

We define 

iAt*B = 


( 88 ) 

(89) 

(90) 

(91) 

(92) 


Note that the phase 0 depends on two quantities: the magnetic fiux normal to the system, and the angle between the Zeeman 
field and the Rashba spin-orbit field. These can both be tuned, and therefore 0 can be viewed as a tunable quantity. 

(«^|(lt;3,effk) ={m\hf,z\n) 

={rn\2iAi*Be-^^-H^-^^pv-^~ne^'P-.p+A^-iSh.v-,^A^AlWxlhi^l+H.c.\n) 


X 7^7p+s7f+s7f + H.c.\n). 


= (m|^ (1 + - 1)55) + 1) + (1 - e*‘^--+^)55+,) 


(93) 
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The eigenstates of interest can be labelled as 

I nz nz nz nz \ 

I ^r—z^r^r-\-x^r-\-x-\-z/ ’ 

Thus, we are interested in the matrix elements: 

We see that the only non-zero matrix elements are those for which {SZ_~y = Si_~, (55)' = —55, (55^^)' 
(55^j_,_^)' = 55^^_,_^. Thus, we need to compute 

(55_,5555+,55+,+,|/i*;3|55_„ -55, -55+„ 55+,+,). 

In terms of these matrix elements, we can then write the effective Hamiltonian: 

hf,3,eff = Y1 (5p-5'5555+,55+,+,|/it.3|55_„ -55, -55+„ 55+,+,) 

{5-} 

^ \qZ qZ qZ qZ \ / I 

\^f_2^f^f-\-x^f-\-x+z/\^f—z’> '^r+x+£ I 

(1 + si55_,) 55 + S2iSy. (1 + 5455+,+,) 


where 


51,52,S3,S4 

3 


5l ,52 ,53 ,S4 = ±1 

= E E c„bcd52_,5^-5£+,5; 

a,d=l,z h,c=x,y 


d 

r-\-x-\-z 5 


^zxxl 




16 

1 

16 


(^zxyl — ^ ^t;3 


52,53,54^^1^3^ 


S 

^ __L J^Si,S2,S3, 

^zyyl — ^ fH-,3 


lSiS2^ 


S 1 S 2 S 3 , 


s 

Clxyz = Yq ^t;3 
s 




18482 


^lyyz 


_ ^ l5i,S2,53,5 

- 10 2^^3 


545253 


^zxyz 


-— /,S1.S2,S3,S4 

■l0Z^^*;3 ®i®4, 

s 

^ ^ ^51,52,53,54 • 

^zyxz —^ iH-,3 ^ 

s 

r 

^zyyz — ^ IH;3 


51 ,52, 53 ,S4-' ^ ^ 

3 ^<51<53'^45 


1 j 8\82 54 5 


51525354, 


C-lxxl — E 

s 

Clxyl ~Yq E 

s 

_ Jl 53,54 • 

fH]3 ^ 

s 

Clyyl “ ~ Ys 55 ^t;£ 


S2,S3,S4 ■ 

3 '*3 


tS2 


S2,S3,S4„ „ 

3 S 2 S 3 


Now, simplifying we get: 


f^sus 2 ,s 3 ,s 4 = ((sg,s2,S3,S4|iIm (e*'^(e - l)(e*‘^-’''+* -f 1)) |si, S2, S3, S4) 

+ S2(si, S2, S3, S4|Re (e*^(e -p l)(e®''’“•'■+*' -|-1)) |si, §2, S3, S4) 

+ 53(51, S2, S3, S4|Re (e®^(e — 1 )(—6*“'’“’''+“’ -|- 1 )) |si, S2, S3, S4) 

+S2S3(si, S2, S3, S4|rtm -P 1)(1 - e*‘'"-'"+»)) |si, §2, S3, S4)) 


Here we have defined 


|si,S2,S3,S4) I SjT — , Si, Sjt S 2 , 5^?+, S 3 , 5^?+,+, S 4 ) 


(94) 

(95) 

-55+,, and 

(96) 


(97) 


(98) 


(99) 


( 100 ) 

(101) 
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The states with the tildes over the s’s indicate that the phase mode has the opposite wave function as compared with the spin 
degree of freedom, as described for the single spin case in Eq. ( [35] ). Let us define 

^31,’s 2 ,33,34 = (si,S 2 ,S 3 ,S 4 |Re +c7i)((T2e*‘^-'"+‘* + 1)) |si, S 2 , S 3 , 54 )- 

^ 3 i;s 2 % 3,34 = (si,S2,S3,S4|«Im(e®®(e"*‘^-'-+(Ti)(cr2e®‘^-’"+* +1)) |si, S 2 , S 3 , S 4 ). (102) 

In terms of ^ we have 


h^U32,3s,34 ^ 


I (■®S1,S2,S3,S4 


+ '®2^s;^S2,33,S4 + '®3^Si^s,,3 


Si ,52,53,54 


We see that we need to compute the following expectation values for the two-spin system: 

^^34,32;± = (55_, = Si, 55 = S2|e±*‘^--|55_, = si, 55 = S 2 ) 

Ws4,32;± = (55_, = Si, 55 = S2|e±'‘"---^ = si, 55 = S 2 ), 

5i,si,s2 = {Sp-z = si,Sp = S2|55_^ = Si, 55 = S 2 ) 

92,34,32 = (55_, = Si, 55 = S2|55:7 ^si, 55 = S 2 ) 

In terms of these expectation values, 


+ S2S3i?^^ ^ 2 , 33 , 34 ) 


(103) 


(104) 


^5 i ,’5''2%3,54 +^15'1,5 i ,52)( c ^2^5354,- +^2,53,54) 

+ -e ^^{Vs^S2,+ + Cri^l,5i,52)(^2'^^S354,+ + ^ 2 , 83 , 54 ) 


1 


-^51,52,53,54 2 ^ (^^5 1 52 , — 3“ 5^1,5 1 ,52 ) (^2 '^84 53 , — “1“ ^2, S 3 , 54 ) 


2^ ('^5152,+ + <^1^1,5i,52)(^2'^S453,+ + ^ 2 , 53 , 54 ) 


(105) 


Note that the eigenstates = si, = § 2 ) can be changed by a phase, which will modify the expressions above. It is useful 
to pick a choice of phase so that, when possible, g^v^w > 0 . 

In Figs. l4][T6 we display some results for the numerical calculation of Cabcd- We see that in all cases, the only appreciable 
couplings are ciabi’, Czbcd and Cabcz are both quite small. This is because the amplitudes have a very weak dependence 

on 5 i and S 4 and thus essentially cancel each other in the sum. 


When Ej = 0, then Fig. 14 shows that at = 0, only cixyi = —ciyxi = 1; for all values of 0, at least two of the four 
couplings cixxi, cixyi, ciyxi, and ciyyi are of the same order. As shown in Figs. T5|[T6 when Ej ^ 0, we can enter the regime 
where only ciyxi is appreciable while all others are much smaller. In fact, we see from the equations above that if the matrix 
element (si, S 2 , S 3 , S 4 |e*'^-’^|si, S 2 , S 3 , S 4 ) is close to one, which is the case for large Ej, then the only appreciable term in 
Cibci will be ciyxi. This is the reason that we need to include Ej ^ 0. The optimal point would be to take Ej to be as large as 
possible; however this would dramatically reduce the Jz coupling, which was calculated in the previous section. Therefore, we 
need to find an optimal point where Ej is non-zero so that only ciyxi is appreciable, while Jz is still large enough. 


2D 


Let us now consider the 2D network shown in Fig. [^. The effective Hamiltonian for this is given by 

H 2 D = ^ ^2s,7 + ^tun, (106) 

I 

where H 2 sg is the Hamiltonian for the Ith unit cell; for each unit cell, this is given by H 2 S above, i^tun contains the horizontal 
tunneling terms, which were analyzed for the case of two horizontally coupled unit cells in the preceding section. For the full 
2D system, it is given by: 


Htun — ^ ^ E tB'4^^y fgpx,f-\-x “ 1 “ H.C.]. 


(107) 
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Cj = l.SCg,Cz = O.SCg,Cz = O.SCg,Ej = 0 


^Ixyl ^Ixxlt ^lyyl 



FIG. 14: Plot of the 16 parameters Cabed, for a, d = l,z and 6, c = x, y, as a function of 0. The only ones that differ appreciably from zero 
are cibci • 


Cj = l.SCg,Cz = O.SCg,Cz = O.SCg,Ej = 0.6E(^j 



FIG. 15: Plot of the 16 parameters Cabed, for a, d = l,z and h^c — y, as a function of 0. The only ones that differ appreciably from zero 
are cibci- We see that for 0 ^ O^tt, ciyxi is much larger than cixyi, cixxi, and ciyyi. 
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Cj = l.SCg,Cz = O.SCg,Cz = O.SCg,Ej = OASE(^j 



FIG. 16: Plot of the 16 parameters Cabcd, for a, c? = l,z and h^c — x^ y, as a function of 0. 
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FIG. 17: Two-dimensional array. 
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Perturbing around the independent unit cell limit = 0, a simple generalization of the analysis of the preceding section 

gives the following effective Hamiltonian, which operates in the subspace spanned by the two effective spin states on each site: 

H2D-Qn = + + ( 108 ) 

R r 

where sums over all unit cells, and R refers to the top spin of each two-spin unit cell. As we have shown in the preceding 
sections, there exist parameter regimes where the additional terms in SH are negligible: 

\\SH\\<^J,,Jy,. ( 109 ) 

H 2 D-,Qff can be recognized to be the Kitaev honeycomb spin model. (531 Specifically, one can perform a tt spin rotation around 
on every other site, which brings the H 2 D]Qff into the form 

R 

which is the more familiar form of the Kitaev model. Here, Jy = —Jx = Jyx- 


Numerical Estimates of Energy Scales 


In Figs. [6p^ we have presented the results of several numerical calculations of the energy spectra of the single effective spin, 
the two spin unit cell, and the couplings constants of the effective spin interaction terms. From Figs. |^-[^ we see that the Jz 
interactions, which couple the vertically separated spins via an interaction S%S% must be on the order of a few percent of the 

R R—z 

Josephson charging energy jCj, in order for the low energy spin manifold to be comfortably separated from the rest of the 
excitations of the system. Figs, pi] show that a finite Josephson coupling Ej is required, so that the horizontal couplings 
will be in the appropriate regime of the Kitaev honeycomb model. 

While we have not performed an exhaustive optimization, our preliminary calculations suggest that the following parameter 
regime is a good one: 


Cj = 1.5C^ Cz = C'z = O.bCg Ej = ^AbEcj ^ = 0, if^l = I^bI ^ 0.1F;gauge. (HI) 


With this choice of parameters, we find that Jz ~ 0.02F^c'j, while the energy cost to the other excited states of the two-spin unit 
cell is approximately ten times as large, F^gauge ~ 0.23F^c'j • This gives a comfortable energy window that separates the low-lying 
effective spin states and the rest of the states of the sytsem. Fig. 16 shows that with this choice of parameters, ciyxi ~ —1.75, 
while cixyi ~ 0.3, and cixxi = ^lyyi = 0, as are all other horizontal coupling terms. This gives almost a factor of approximately 
6 between the horizontal couplings that we want and the undesired ones. In terms of absolute energy scales, we have: 


141 = 141 = 1.75|fAiB|/4 


gauge- 


( 112 ) 


If we set |f AI = = 0.2F^gauge to ensure the single electron tunneling processes are suppressed relative to the second order 

process, we find \Jx\ = \Jy \ = 1.75 x 0.04F^gauge = O.OlfiF^Cj. which is almost the same order as the Jz estimate above. 

Therefore, we see that to get an effective spin model whose dominant interactions are the Kitaev interactions, while all other 
interactions are suppressed, we can get energy scales that are roughly in the range of a few percent of the Josephson charging 
energies Ecj . To get a large energy scale, then, we wish to use a physical setup with the largest possible Josephson charging 
energy e^/Cj, which can also simultaneously accommodate a Josephson coupling Ej « O.be^jCj. 

Typical Al-Al^^Oi-a^ Josephson junctions have Josephson charging energies on the order of Ecj ~ 1 K, which can therefore 
give interaction strengths Jz^Jx^Jy ~ 20 mK. 

Josephson junctions made from gated semiconductor wires, such as Al-InAs-Al junctions, can yield much larger Josephson 
charging energies, because the distance between the superconductors (ie the length of the nanowire junction) can be much larger. 
For example, let us consider InAs wires with radius r and a distance d between the A1 superconductors. For r = 20 — 60 nm 
and d = 100 — 450 nm, critical supercurrents Jc = 1 — 135 nA have been measured, El which corresponds to Ej = hld'le « 
0.05 — 3 K. If we consider d = 100 nm and the superconductor consisting of a wire of A1 epitaxially grown on the InAs nanowire 
with total radius 100 nm, we can estimate Cj = ereo7r(100nm)^/(100nm). Taking = 15 for InAs, this implies a charging 
energy d jCj « 40 K. This can be further reduced by increasing the radius or decreasing d. 

Interestingly, devices with d = 30 nm have also been fabricated, and have been reported to yield critical currents as high as 
Ic = 800 nA, for InAs nanowires with radius r = 40 nm. |[64l This corresponds to Ej ^ 40 K. If we assume the parallel plate 
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capacitor formula for a radius 40 nm and d = 30 nm, we would get jCj ^ 75 K. However, it is not clear whether such a high 
supercurrent is due to unwanted parasitic effects that are introduced during the fabrication process. 

To put this on a somewhat more theoretical footing, consider that the supercurrent is typically given by 

/ci^AT = 7rA/2e, (113) 

where i?Ar is the normal-state resistance of the junction, and A is the superconducting gap. For a semiconducting wire with Nc 
channels, this implies 

p- 2 

Ej = hlj2e = ^ 7 ATVex = ^^c/8 (114) 

4 ri 

where the conductance is e^//i per channel. For Al, with A = 1.2 K, this implies that Ej = 0.15Ac K 

The above considerations, and in particular the experimental measurements, suggest that it could be possible, with Al-InAs-Al 
junctions, to get to a regime where Ecj = 2Ej ^ 5 — 10 K. This would then imply 

J„J^,J^^0.1-0.2K. (115) 

Note that Nb is also a candidate material that can be used in these setups, instead of Al. Indeed, Nb-InAs-Nb Josephson junc¬ 
tions have been fabricated and measured. (651 While the use of Nb presents certain technical obstacles for fabricating the required 
semiconductor-superconductor heterostructures, it has the advantage that the superconducting gap is much larger, ANb ~ 9 K. 
This implies that the energy scales considered above will be a factor of Aa\ = 7.5 larger if Nb is used instead and good 
contact can be made between the Nb and the InAs. So far such an enhancement in the critical supercurrent has not been observed 
due to contact quality, but there are no fundamental obstacles to improving this contact quality and thus achieving this factor of 
7 — 8 enhancement. 

This would then suggest the theoretical possibility 

J„J^,J^^0.5-2K. (116) 

Superconductors with even larger gaps, such as NbTiN, could potentially yield even larger energy scales. 


Quantum Phase Slip Based Implementation 
Single spin 


We now consider an alternative possible architecture, which utilizes the physics of quantum phase slips to engineer an effective 
Kitaev spin model. The main building block of this architecture is a set of four Majorana fermion zero modes, as shown in Fig. 
[T^. Each pair of Majorana zero modes arises from the endpoints of a spin-orbit coupled semiconductor nanowire, proximity 
coupled to a superconducting island. Each pair of Majorana fermion zero modes gives rise to two degenerate states, associated 
with whether the fermion parity on the island is even or odd. The four Majorana zero modes will be labelled as 7 ^, 7 ^, 7^, and 
7 ^, as shown in Eig. 


18 \. 


The effective Hamiltonian for the system shown in Eig. is 


Ass = ^ ABdG[Ao 

j=A,B 


- Ej cos{ipA - (Pb) 


cos(^i) + J ^ 


Ej^ 2_^ cos 
i=A,B 


Q^c^^QJ 


(117) 


'd=A,B 


Here, cpj for j = A^B is the superconducting phase on the A and B islands, Is the BdG Hamiltonian for 

the nano wire on the jth island, where | Aq^ | is the proximity-induced superconducting gap on the jth nano wire at zero magnetic 
field. Qj is the excess charge on the jth superconducting island - nanowire combination; it can be written as: 

Qj = e{-2id^. -h Nj - rioffj), (118) 

where represents the number of Cooper pairs on the jth superconducting island, Nj = f is the total number of 

electrons on the jth nanowire, and Uoffj is the remaining offset charge on the jth island, which can be tuned continuously with 
the gate voltage Vg . The capacitance matrix is 


(Cg Cj + Cr —Cj \ 

V -Cj CgECjE Cr) 


(119) 
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FIG. 18: (a) Single effective spin. The A and B islands are each coupled to a large superconductor with Josephson coupling Jr. This is to 
emulate the embedding of the single spin into a larger network, (b) ID chain. Ia and indicate electron tunneling, (c) A single plaquette of 
the 2D network. The Josephson coupling between different spins is Ej ^, while Ej-^ is the Josephson coupling between the A and B island of 
a single spin. 


As in Sec. | we perform a unitary transformation U = to decouple the phase cpj from the 

fermions 2 pj in iTBdG- Here, umj = 0,1 is the occupation number of the pair of Majorana zero modes on wire j. It is given in 
terms of the Majorana zero modes as 

nMA = (1 + i7''7*)/2, 

nMB = (1 + h''1^)12. (120) 

Under this transformation, the charge Qj transforms as: 

Q'j = U^QjU = e{-2id^^ 

Thus, taking Hgs U^HssU, we obtain 

=U^H,,U = Y, ffBdG[Aoi, V’], - Ej COsi^A - ^b) - Ej^ Y ^ S Q'iC^^Q'j 

3 i=A,B i,j=A,B 

= Y EBdG[Aoj,Aj,ipj] - Ej cos((^_) - 2Ej^ cos{ip+/2) cos((^_/2) 

3 

+ \iQ+fic^\ + cJb) + - Ca^), ( 122 ) 

where we have used Caa = Cbb, and we have set 

Q±=Qa±Qb. (123) 

In what follows we will drop the prime superscripts for convenience. 

It is helpful to consider the Lagrangian for this system, which is given by 

Ev = ^ - nos)i^i + Ej cos{i^a - ^b) + Ej^ Y cos{i^i) 

i=A,B 


(124) 
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We consider the limit where 


Ej,Ej^->e‘^Cj}, (125) 

in which case (fA, are pinned, while the conjugate variables Na, Nb are highly fluctuating. For energy scales below 
Ej^ Ej^, the effective Hamiltonian of this two island system takes the form 

Hob = + Cb*7 ^7‘ - CabAi^iA* (126) 

The first two terms are due to quantum phase slip events where either cpA or cps change by 27r. The last term is due to the 
quantum phase slip event where both islands A and B collectively change their phase by 27r, relative to the phase of the 
reservoir. These phase slip processes effectively measure the fermion parity of the region undergoing the phase slip, which can 
be expressed in terms of the Majorana fermion modes. 

The quantum phase slip amplitudes, CaXb, Cab^ can be computed using the standard instanton calculation, in the limit of 
dilute instantons. To leading order in ^ we estimate this to be 

Cb = -cos{,n.,B)8£o„,/? 

V TT y 2E(j^^ J 

Cab = - cos(7rnoff+)8i;c^^ ' e” (127) 

V TT y J 

where Ec,^ = 5^, and Ec^.^ = t(<^aa + C'ab)- 

We assume that the direct coupling between these different Majorana zero modes, which is generated by electron tunneling 
between the two ends of the wires, is much smaller than all other energy scales in the problem, and can therefore be ignored. 
Observe now that if we set 


noff7 = l/2, / = A,B, (128) 

then the effective Hamiltonian is simply 

H,«,ss = -UbiAAA\ (129) 

which enforces the constraint 

for states with energies much less than (ab- The ground state subspace of this system is therefore doubly degenerate and acts 
like a spin, with 

(131) 

Slightly tuning the offset nofM and/or UoffB away from 1/2 acts like a Zeeman held hzS^. 


ID Chain 

Next, let us consider a ID chain of the ^4, 5 superconducting island pairs introduced in the previous subsection. The coupling 
between the effective sites of the chain consists of the Majorana fermion tunneling terms shown in Fig. [T^, with tunneling 
amplitudes Ia and Ib- In addition to the single particle tunneling terms between the Majorana zero modes, there are pair 
tunneling terms that induce Josephson couplings Ja, Jb between the islands. Finally, there are cross-capacitances between the 
different superconducting islands; we will consider the regime where the charging energies due to these cross-capacitances are 
much smaller than the tunneling and Josephson couplings, and can therefore be ignored. The effective Hamiltonian of the chain 
is therefore 


Hid — ^ ^ [^r ~l~ 


(132) 
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where each hf is given by the Hamiltonian Hss described in the previous section, and 

hf^f-\-x ~ ^ “ 1 “ tB'^yfM^xf-l-x 

Ja iffj^x^A) Jb COs((/9-^^^ ^f+x^B^ (133) 

Note that the pair tunnelings Ja and Jb in this setup are generated by pair tunneling between the Majorana zero modes, and 
thus Ja < t a, Jb < ^b- Performing the unitary transformation U in the previous section, and then setting 

= (134) 


as in the previous section, we get: 

= [t^e“*(l + ^P>^>r)¥’Ar/2 + i(l-Fp,A,r+£)^Ar+x/2^t^t^^ 

+ + H.C.] 

Ja ^r-\-x,A) Jb COs((^f?^B ^f+x,B\ (135) 

Now, since we are in the limit of large Josephson couplings, all of the phases of the superconducting islands can be set equal to 
each other, which we can set to zero without loss of generality: 

^Ar = ^Br=0, (136) 

with corrections coming from instanton events. Thus, becomes 

hr,r+x = [2iIm(B)7A7s + 2«Im(iB)77?+s]- (137) 

We wish to consider the limit where 


tA,tB < Cab- 


(138) 


In this case, the single tunneling events are suppressed; perturbing to second order in we obtain an effective Hamiltonian: 


r,r-\-x 


4Im(fA)Im(tB) V t X t 

— A hhf+sA+x 


Therefore, the effective Hamiltonian of the ID chain is given by: 


(139) 


Hid = '^[CAi'lhe + Csilhf- “ UBlhhflf “ 


(140) 


If we assume for simplicity that Ja = Jb = J, then the quantum phase slip amplitudes are essentially as given in Eqn. |147[ 
with Ej^ = 2 J, and Cr being twice the capacitance across the semiconductor wires. 


2D Network 

Now we would like to assemble the ID chains described above into a two-dimensional network. We consider the network 
shown in Fig. [T^, which effectively forms a brick (honeycomb) lattice. Each unit cell of the lattice consists of two pairs of 
superconducting islands: two A islands and two B islands, separated vertically from each other. The Hamiltonian for the system 
can be written as 


H 2 D=E^ 2 e,R + ^t<^«’ ( 141 ) 

R 

where R is the location of the top spin of each unit cell, is thus a sum over unit cells, and the Hamiltonian for each unit cell 


QiCjjQj Ej^COs{ipj^^ ^Br) J^BR-z) 

1=1 ^ I,J=1 

— Ejz ^^^{^BR ~ ^AR-z)' 


(142) 
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The 4 X 4 capacitance matrix now is 


/Cg+Cj -Cj 0 0 \ 

^ _ —Cj Cg + Cj + Cj^ ~Cjz 0 

0 -Cj, C, + Cj + Cj, -Cj 

VO 0 -Cj Cg + Cj) 

The tunneling Hamiltonian i^tun is: 

Ht ^ H.C.] 


(143) 


(144) 


We have ignored the capacitance between horizontally separated islands, as they are assumed to be far enough apart that their 
capacitance is negligible. 

Based on the analysis of the single spin case and the ID chain, we can now immediately see that the low energy effective 
Hamiltnian for this system can be written as 


^2D;eff = + i^2 + ^t;eff, 

where Hi consists of single island phase slips: 


R 


H 2 consists of two-island phase slips: 

H, =-Y. if" zhhH - 

r R 

- Cx'^lhhhx'rp+x - C|^ 7 p 7 f 7 f+s 7 f+* 


(145) 


(146) 


(147) 


These phase slip amplitudes are given in terms of the Josephson couplings and charging energies: 


C 4 = - y I ( e- 

cl = - C«s(™.„,g)8£c.. f- (^1^) 3-)/^o.. 


Cf . = - cos(7rn ffo s_ J8^C44\/ z ( e V^(Ej,)/Ec, 


2 / E 


TT V 2K 


044 


3/4 


(148) 


aAB _ _ 
Sr 


cr 


where we have defined 


COs(7r(77/ofM,r H“ '^offB,r))8-£/C 
cos(7r(noffB.f + no{{A,i^-z))8Ec 

= 


2 f Ej 


3/4 


J^] -^S(Ej^)/Ec^ 

^^^Ix\2Eca 


-^AB 

2 ( 2Ej, 

^^^'ttUEca 


3/4 


- j8(2Ejj)/Ec£ 


2 2C„ + Cz 


(149) 


^Cba = 


2 2Cg + 2Cj 


(150) 
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The phase slip amplitudes Cx^ approximately given by products of the single island phase slips, 

/•AA ^ /-A /-A 
/-BB ^ /-B/-B 


(151) 


because the horizontal Josephson couplings are negligible. 

All other two-island phase slips, and collective phase slips of more than two islands, can be ignored, as their amplitudes are 
exponentially suppressed relative to the terms considered here. Finally, i^tun includes the Majorana tunneling terms Ia and Ir'- 


f 


(152) 


If we set nofM,r = =1/2, then the single island phase slips vanish. The effective Hamiltonian becomes 


H2d-,^s = - y] Cp^ihhpip 


- Y' MB X y z^ 

IrIrIr-z^R-z 


4Im(f^)Im(fB) ^ _ 

Ir Ir ir-\-x ir-\-x 


Cab 


If we further consider the regime where 

0 ^ » Cz^,2Im(f^),2Im(fB), 

we see that the system can be described by an effective spin model: 

ff2D;eff = ^ E ^Phx, 


(153) 


(154) 


(155) 


with 


R 


T — 


Jyx 


4Im(/^)Im(/B) 

Cab 


(156) 


Upon rotating every other spin around the .2 axis by 7r/2, the above Hamiltonian can be put into the more familiar Kitaev 
form: 


^f2D;eff= E E E 

{ij) =z-link ( ij) =x-link ( ij) =y-link 


(157) 


Realizing the Ising topological order 

As we have shown, the physical architectures described above can give rise to an effective realization of the 2D Kitaev 
honeycomb spin model: 


77k = J. E +JyY. (158) 

x-links y-links z-links 

It is well-known that this model is proximate to a non-Abelian topological state with Ising topological order. There are a number 
of known ways to access the non-Abelian state. One way to access the non-Abelian state is to apply a small effective Zeeman 
field: 


SH = }2 hxS^ + hySl + hzSZ 


(159) 


As explained in Sec. [] The hz term above can be generated by tuning the gate voltage on each of the superconducting islands. 
The hx and hy terms are more difficult, but possible, to generate as well. They require connecting the Majorana zero modes with 
additional semiconductor wires, as shown in Figs. [^, to allow for electron tunneling as shown. 

A second way of realizing the Ising topological order is to consider effectively the same model, but on a different lattice (see 


Fig. 19), as proposed by Yao and Kivelson. As a spin model, the Kitaev Hamiltonian, Eq. (158) is time-reversal invariant. On 
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FIG. 19: Micro-Architecture for the Yao-Kivelson realization of the model. 


the lattice structure proposed by Yao-Kivelson the ground state spontaneously breaks time reversal symmetry, yielding a ground 
state with Ising (or its time-reversed partner, Ising) topological order. In the completely isotropic limit where all couplings are 
equal to J, the energy gap of the Ising state is also equal to J. Interestingly, disorder in the spin couplings can actually be 
beneficial and can enlarge the region of stability of the Ising phase. (661 

By adding a small effective time-reversal symmetry breaking perturbation to the spin model, we can tune whether the topolog¬ 
ical order is Ising or Ising, and avoid having domains of either, as would be realistically expected in the case where the effective 
time-reversal symmetry is broken spontaneously. A Zeeman term by itself is insufficient. We can consider either a 

Zeeman term that includes both hz and hy. Or, in order to avoid requiring hy or hx terms, which are more difficult to generate, 
we can make use of the smaller perturbations that are naturally generated in the first charging energy based 

implementation presented above. 


Genons 


Ising X Ising topological order 


Given a microscopic architecture to realize a quantum state with Ising topological order, one can then consider designing two 
independent copies of such a state (referred to as the Ising x Ising state) by utilizing present-day nanofabrication technology to 
create short overpasses among different superconducting wires, as shown in Fig. 


Creating genons 


A genon in an Ising x Ising state can then be realized by modifying the overpass connections to create a segment along which 
the connections among the horizontal chains is twisted, as shown in Fig. These segments effectively create branch lines that 
connect one layer to the other, and vice versa. The end-points of the segments realize exotic non-Abelian twist defects, which 
have been referred to as genons. The topological degeneracy of the system in the presence of the genons mimics that of a single 
copy Ising system on a surface of non-trivial topology. (26l[3ll|5ll 

Technically, there are three topologically distinct types of genons, which we label as Xj, and X^. lfSSl Physically, they 
correspond to whether a I, t/;, cr particle, from either layer, is bound to the genon. The genons have the following fusion rules: 
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FIG. 20: Two essentially decoupled copies of the capacitance-based model can be created by using short overpasses to couple next nearest 
neighbor chains. 


Xix{I,i;) = Xix{ij,I)=X^ 
Xi X (I, a) = Xix (cr, I) = Xc, 


(160) 
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-# — #■ 

(b) 



FIG. 21: Architecture for creating genons in the effective spin model. The location of the genons corresponds to the end point of the branch 
lines, and are marked by green circles, (a) Full lattice dislocation, (b) The vertical bonds that skip two chains along the branch cut are removed. 
The fact that these vertical bonds can be removed follows the analysis of Ref. |67l which showed in a different context that half of the branch 
cut is sufficient. 


It follows from the above that Xj, have quantum dimension 2, while X^j has quantum dimension 2\/2. 


Effective braiding of genons and the topological tt/S phase gate 

In Ref. [3TJ it was shown that the braiding of genons can be used to realize a topologically protected tt/S phase gate. In 
the present system, the braiding of genons is complicated by the fact that it is difficult to continuously modify the physical 
location of the genons to execute a braid loop in real space. Fortunately, this is not necessary, as the braiding of the genons can 
be implemented through a measurement-based approach. To do this, we require that it be possible to measure the joint fusion 
channel of any pair of genons and project it into either the (I, I) channel or the (t/;, 2 p) channel. 


Measurement-based braiding of genons 


Importantly, the braiding of the genons can be achieved without moving them continuously around each other in space, but 
rather through tuning the effective interactions between them. Specifically, what is required to braid two genons is the ability to 
project the fusion channel of pairs of genons onto an Abelian charge sector. 

In order to implement the tt/S phase gate, we wish to start with two pairs of genons, labelled 1,4, and have the ability to 
braid genons 2 and 3. In order to do this, we use an ancillary pair of genons, labelled 5 and 6 . The braiding process is then 
established by projecting the genons 5 and 6 onto the fusion channel then the genons 5 and 3 onto the fusion channel 635 , 
the genons 5 and 2 onto the fusion channel 625 , and finally again the genons 5 and 6 onto the fusion channel 655 . 

We will asssume here that b^e = b'^Q. If the genons 5 and 6 are created out of the vacuum, then it will in fact be natural to 
have 656 = b'^Q = (I, I). In this situation, we can derive the resulting braid matrix for the genons, following the results of Ref. 
[ 59 I When 625 = 635 , the braid matrix for genons 2 and 3 is given by 


R23 = e 



0 

0 

gi7r/8 


1 0 

0 -1 
0 0 


(161) 
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FIG. 22: Schematic for measurement based braiding of genons. In order to effectively braid genons 2 and 3, the following protocol is 
performed, (a) The fusion channel of genons 5 and 6 is projected onto the any on b^e. (b) The fusion channel of genons 5 and 3 is projected 
onto 653 . (c) The fusion channel of genons 5 and 2 is projected onto b 52 - (d) The fusion channel of genons 5 and 6 is projected onto b'^Q. 


where is an undetermined, non-topological phase. In other words, the state obtains a phase of ±1 or depending on 

whether the fusion channel of genons 2 and 3 is (I, I), (t/^, t/;), or (a, a). If instead 625 7 ^ ^ 35 , then 

/-10 0 \ 

i?23 = 0 1 0 (162) 

\ 0 0 


Physical implementation of projection of pairs of genons onto (I, I) or (^, fusion channels 

When two genons are separated by a distance L, the effective Hamiltonian in the degenerate subspace spanned by the genons 
obtains non-local Wilson loop operators: 

-^genon ^('0,'0) (^) 3” ^(cr,cr) ^^(cr,cr) (^) H.C. (163) 

describes the exchange of a particle between the two genons, which equivalently corresponds to a (I, t/^) 

or particle encircling the pair of genons. Similarly, describes the exchange of a (cr, a) particle between the 

two genons, which equivalently corresponds to a (I, cr) or (cr, I) particle encircling the pair of genons. oc and 

oc are the tunneling amplitudes, with and Ccr being the energy gaps for the t/; and cr particles, and some 

appropriate velocity scales. 

Next, let us suppose that the pair of genons shown in the figure fuse to the quasiparticle ( 6 , b). The outcome of the process 
where a ( 1 , a) quasiparticle encircles a topological charge ( 6 , b) is determined by the topological S matrix of the Ising phase, 
and is given by Sah /where 


/ 1/2 1 / V 2 1/2 \ 

S=\l/^/2 0 -l/\/2 , (164) 

V 1/2 -1/V2 1/2 J 

and where the entries are ordered I, cr, t/;. In other words, the eigenvalues of W(^a,a) {C) are given by Sah / Shi, where ( 6 , b) is the 
fusion channel of the two genons connected by the path C. 

The ground state of i^genon, which depends on therefore corresponds to a definite fusion channel for the pair of genons 
involved. We can distinguish the following possibilities: 

1 . I ter I > 1^1 and to- > 0. The ground state subspace of i^genon corresponds to the case where the two genons have fused to 
the (t/;, t/;) channel. 
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FIG. 23: Ising xlsing system in the presence of two holes (three gapped boundaries). This is effectively equivalent to the Ising x Ising 
system with 6 genons. The topologically robust operation described in Fig. [^can be adapted to this case, by projecting the topological charge 
through the loops shown to be equal to 5i,..., ^ 4 . The solid line indicates that it is in the ’’top” layer, while the dashed line indicates that it is 
in the ’’bottom” layer. 


2 . I ter I > |t^| and to- < 0 . The ground state subspace of i^genon corresponds to the case where the two genons have fused to 
the (I, I) channel. 

3. I ter I < |t^|, t^ < 0, to- > 0. The ground state subspace of i^genon corresponds to the case where the two genons have 
fused to the (t/;, channel. 

4. I ter I < |t^|, t^ < 0, ter < 0. The ground state subspace of i^genon corresponds to the case where the two genons have 
fused to the (I, I) channel. 

5. I ter I < |t-01, t^ >0. The ground state subspace of i^genon corresponds to the case where the two genons have fused to the 
(cr, cr) fusion channel. 

We see that there is only one possibility that needs to be prevented, which is the last one, where \ta\ < \t^\, > 0. To do 

this, we flip the sign of by flipping the sign of the coupling along a single link of the shortest path that connects the genons. 
Moreover, note that one can also pick the precise path C along which the quasiparticles tunnel by depressing the gap along that 
path, which can be done by decreasing the couplings along that path. 


Ising xlsing topological order and gapped boundaries 


In the above, we have suggested creating genons in the Ising x Ising topological state. The braiding of the genons, which can 
be performed in a measurement-based fashion, can be used for the topologically protected 7 r /8 gate. Here we note that the same 
topologically robust transformations can also be achieved with the Ising xlsing topological state in the presence of multiple 
disconnected gapped boundaries, where Ising refers to the time-reversed conjugate of the Ising state. 

The main observation is that the Ising x Ising state in the presence of n disconnected gapped boundaries can be viewed as 
effectively a flattened version of a single Ising state on a genus g = n — 1 surface. This is similar to the fact that an Ising x 
Ising state, in the presence of n pairs of genons, can also be effectively mapped onto a single Ising state on a genus g = n — 1 
surface. The measurement-based braiding protocol for the genons of the Ising x Ising state can be readily adapted to the case of 
the Ising xlsing with gapped boundaries. Below we will briefly sketch this adapted protocol, leaving a more detailed discussion 
for future work. 

Let us consider the Ising xlsing state in the presence of two disconnected gapped boundaries, i.e. on an annulus. This is 
equivalent to a single Ising state on a torus, similar to the case of the Ising x Ising state with four genons. In order to carry out 
an effective Dehn twist in this effective torus, we use an ancillary gapped boundary, as shown in Fig. This system is now 


equivalent to a genus 2 surface. In order to effectively carry out the Dehn twist, we perform a series of projections along various 
loops in the system. We consider the loops shown in Fig. [23]\-D, and sequentially project the topological charge through those 
loops to be 61 , • • • , 64 . This is precisely analogous to the genon braiding case described in Fig. where the fusion channel of 
genons was projected onto 656 , 635 , 625 , and 655 . Thus if we take bi = 64 and 61 , • • • , 64 to be all Abelian topological charges, 
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the equivalence between the Ising x Ising system with genons and the Ising x Ising with gapped boundaries implies that we will 
have effectively carried out the desired operation. 

As in the case of the genons in the Ising x Ising state described in the previous sections, these projections can effectively be 
implemented by reducing the gap for quasiparticle tunneling alon g the various loops as required. _ 



